AN EFFICIENT METHOD FOR MULTIOBJECTIVE OPTIMAL CONTROL AND OPTIMAL CONTROL SUBJECT TO INTEGRAL CONSTRAINTS

We introduce a new and efficient numerical method for multicriterion optimal control and single criterion optimal control under integral constraints. The approach is based on extending the state space to include information on a ＂budget＂ remaining to satisfy each constraint; the augmented Hamilton-Jacobi-Bellman PDE is then solved numerically. The efficiency of our approach hinges on the causality in that PDE, i.e., the monotonicity of characteristic curves in one of the newly added dimensions. A semi-Lagrangian ＂marching＂ method is used to approximate the discontinuous viscosity solution efficiently. We compare this to a recently introduced ＂weighted sum＂ based algorithm for the same problem [25]. We illustrate our method using examples from flight path planning and robotic navigation in the presence of friendly and adversarial observers.     

Robin型离散Schwarz波形松弛算法的收敛性分析

Schwarz波形松弛(Schwarz waveform relaxation,SWR)是一种新型区域分解算法,是当今并行计算研究领域的焦点之一,但针对该算法的收敛性分析基本上都停留在时空连续层面.从实际计算角度看,分析离散SWR算法的收敛性更重要.本文考虑SWR研究领域中非常流行的Robin型人工边界条件,分析时空离散参数t和x、模型参数等因素对算法收敛速度的影响.Robin型人工边界条件中含有一个自由参数p,可以用来优化算法的收敛速度,但最优参数的选取却需要求解一个非常复杂的极小-极大问题.本文对该极小-极大问题进行深入分析,给出最优参数的计算方法.本文给出的数值实验结果表明所获最优参数具有以下优点:(1)相比连续情形下所获最优参数,利用离散情形下获得的参数可以进一步提高Robin型SWR算法在实际计算中的收敛速度,当固定t或x而令另一个趋于零时,利用离散情形下所获参数可以使算法的收敛速度具有鲁棒性(即收敛速度不随离散参数的减小而持续变慢).(2)相比连续情形下所获收敛速度估计,离散情形下获得的收敛速度估计可以更加准确地预测算法的实际收敛速度.     

Nonlinear PDE Model for Image Restoration Using Second-Order Hyperbolic Equations

In this article we consider a novel nonlinear PDE-based image denoising technique. The proposed restoration model uses second-order hyperbolic diffusion equations. It represents an improved nonlinear version of a linear hyperbolic PDE model developed recently by the author, providing more effective noise removal results while preserving the edges and other image features. A rigorous mathematical investigation is performed on this new differential model and its well-posedness is treated. Next, a consistent finite-difference numerical approximation scheme is proposed for this nonlinear diffusion-based approach. Our successful image denoising experiments and method comparisons are also described.     

On the proximinality of compact operators over spaces of measurable functions

The main results in this paper are, the characterization of all the σfinite positive measures μ and v, for which K(L₁(μ, L₁(v)) proximinal in K(L₁(μ, L₁(v)) all the a-finite positive measures fi and v, for which K(L_∞(μ), L_∞(v)) is proximinal in L(L_∞μ)), L_v(v)) and all the compact Hausdorff spaces Q, for which K(C(Q), L_∞(μ)) is proximinal in L(C(Q), :_∞(μ))     

A blow-up theorem for a discrete semilinear wave equation

In this paper, a discretization of a semilinear wave equation whose nonlinear term is a power function is investigated. It is shown that, when a condition on the initial value problem, similar to that governing the existence of blow-up solutions for the original continuous equation is met, the newly introduced difference equation has blow-up solutions with characteristics corresponding to those of the blow-up solutions for the original equation.     

Tropical discretization: ultradiscrete Fisher–KPP equation and ultradiscrete Allen–Cahn equation

A systematic approach to the construction of ultradiscrete analogues for differential systems is presented. This method is tailored to first-order differential equations and reaction–diffusion systems. The discretizing method is applied to Fisher–KPP equation and Allen–Cahn equation. Stationary solutions, travelling wave solutions and entire solutions of the resulting ultradiscrete systems are constructed.     

高阶控制增稳系统等效拟配参数时域在线辨识方法

现代控制增稳飞机是高达几十阶的复杂系统,将其降阶为特定形式的低阶等效系统才能与已有飞行品质规范进行比较,进而预测和评价飞机飞行品质.从时域角度出发,提出了一种高阶控制增稳系统等效拟配参数在线辨识方法.分析了时域等效系统拟配的基本方法,提出了等效系统实时在线拟配的系统框架思路,基于Tustin变换法推导了由高阶系统时域响应获取四阶等效系统参数的数学公式,应用最小二乘法在MATLAB中实现了算法编程,并用试飞数据验证了所提方法的可行性.计算结果表明,所提方法可以快速准确地实现高阶控制增稳系统的等效系统参数辨识.     

A comparison of spatial discretization schemes for differential solution methods of the radiative transfer equation

A comparison of discretization schemes required to evaluate the radiation intensity at the cell faces of a control volume in differential solution methods of the radiative transfer equation is presented. Several schemes developed using the normalized variable diagram and the total variation diminishing formalisms are compared along with essentially non-oscillatory schemes and genuinely multidimensional schemes. The calculations were carried out using the discrete ordinates method, but the analysis is equally valid for the finite-volume method. It is shown that the S schemes of the genuinely multidimensional family perform quite well, particularly in problems with discontinuous radiation intensity fields. However, they are time consuming, and so they do not always become more attractive regarding the trade-off between accuracy and computational requirements, in comparison with other high-order schemes that, although being less accurate, are also more economical.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.